Nonlinear conductance of long quantum wires at a conductance plateau transition: 

Where does the voltage drop? 



T. Micklitz\ A. Levchenko^, and A. Rosch^ 
^ Dahlem Center for Complex Quantum Systems and Institut fiir Theoretische Physik, 
Freie Universitdt Berlin, 14195 Berlin, Cermany 
^Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 
"^Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77, 50937 Cologne, Germany 

(Dated: February 19, 2012) 

We calculate the linear and nonlinear conductance of spinless fermions in clean, long quantum 
wires where short-ranged interactions lead locally to equilibration. Close to the quantum phase 
transition where the conductance jumps from zero to one conductance quantum, the conductance 
obtains an universal form governed by the ratios of temperature, bias voltage and gate voltage. 
Asymptotic analytic results are compared to solutions of a Boltzmann equation which includes the 
effects of three-particle scattering. Surprisingly, we find that for long wires the voltage predominantly 
drops close to one end of the quantum wire due to a thermoelectric effect. 

PACS numbers: 71.10.Pm, 72.10.-d, 72.15.Lh 



Introduction. — A clean quantum wire with adiabatic 
contacts is characterized by a quantized conductance, 



G 



Gn with Gn 



' /h. The integer n describes the 



number of conduction channels (including spin). The 
conductance quantization is closely related to charge 
quantization and survives (for sufficiently low tempera- 
tures T) even in the presence of interactions [iHl] as long 
as momentum relaxation by Umklapp scattering can be 
neglected [1, 

The transition from one conductance plateau to the 
next is an example of a quantum phase transition with- 
out order parameter, where only a topological property, 
the number of conducting channels, changes. While the 
thermodynamics of this quantum critical point (QCP) 
is quite well understood [o-'fl'], a theory of the quantum 
critical conductance is much more challenging: How does 
the conductance change from one conductance plateau 
to the next at low but finite T? We will answer this 
question both in the linear and non-linear regime for the 
most simple situation, i.e., the transition from n — Q 
to n = 1 for spinless fermions with finite-ranged interac- 
tions. Here the QCP describes the transition from zero to 
a finite fermion density. As interactions are irrelevant at 
this QCP, thermodynamic pro per ties are well described 
by non-interaction fermions [sl llOf but transport in long 
wires is still governed by collisions. Relaxation by coUis- 
sions and non-equilibrium dynamics in one-dimensional 
(ID) systems have recently moved into the focus of the- 
oretical [ll-18| and experimental 1^21| research. 

A general question is where the voltage drops when a 
finite current is driven through a clean ID quantum wire 
by applying a bias voltage V . For diffusive (multichan- 
nel) quantum wires one expects a linear drop of the volt- 
age (i.e. of the electrochemical potential) across the wire 
while for non-interacting, ballistic quantum wires the 
voltage drop occurs only close to the two contacts [2^ . 
In clean interacting quantum wires with low fermion den- 



sity (and therefore negligible Umklapp scattering) the dc 
conductivity is infinite for an infinitely long wire due to 
momentum conservation, a{T) = oo. The vanishing re- 
sistivity strongly suggests that there is again no voltage 
drop inside the wire. 

A recent series of papers 

flSIHMliS, which studied 
the role of equilibration in long but finite quantum wires 
of length L, found that in the linear response regime, 
1^ — )■ 0, there is a linear drop of voltage [3| along the 
wire. We resolve this apparent contradiction to (t(T) = 
oo by noting that the limits V ^ Q and L — )■ oo do 
not commute. The drop of voltage is governed by a new 
length scale which diverges for — )■ 0. For L <C a 
linear drop of voltage occurs. In the opposite limit, L 3> 
Iv, however, the voltage drops only within a distance of 
Iv of the contacts. Surprisingly, the voltage drop is not 
symmetrical and occurs predominantly only at one of the 
two contacts! 

Previous work on equilibrated quantum wires P, 0, 



24 l26l . |27| focused on the limit T ep, where Ci;- is 
the Fermi energy. As scattering processes equilibrating 
left-moving and right-moving fermions involve the bot- 
tom of the band [2|, [18|, [ii-lii], they are exponentially 
suppressed and therefore the corresponding equilibration 



length is exponentially large. 



e^l'^ , where A e^? 



for weak interactions (thermal equilibration relevant for 
heat conductance occurs on shorter length scales 29]). 
For L 3> ^oq it was found that the quantized (linear) 
conductance obtains corrections of order [T/ /S.)^. Large 
effects can therefore be expected close to the conductance 
plateau transition, where A ~ T, as studied in this pa- 
per. 

Model. — We consider ID spin-polarized electrons with 
quadratic dispersion Cp — interacting via a short 
range potential. Close to the QCP, where filling of the 
first subband becomes small, interactions are strongly 
irrelevant in the renormalizati;,n group sense iH and a 



2 



single electron description becomes approximately valid. 
To study equilibration and its effect on transport we may 
thus use the Boltzmann equation 



Vpdxf: 



(1) 



where fx,p is the quasiclassical distribution function, 
p/m the velocity, and the collision integral 1"^°^ 



Vp 

describes collisions. The contacts of the quantum wire 
to the leads at a; = ±-| induce boundary conditions for 
electrons moving into the quantum wire 



1 



,p>o 



1 



fx- 



1 



f,P<0 



1 



(2) 



where S^J^ — Sp — ^J■T eV/2 with fi — at the QCP. Here 
we assume adiabatic and ballistic contacts, i.e., contacts 
which are smooth compared to the electronic wavelength 
but short in comparison to the scattering length. 

In ID systems energy and momentum conservation 
severely restrict the phase space available for scattering: 
in a two-particle process, two particles of equal mass can 
only exchange their momenta d, [l^] which leaves fp un- 
changed. One therefore has to study the effects of three 
particle collisions [H, 28| described by 



E W^123'''[/l/2/3(l-/l')(l-/2')(l-/3') 

i 

-/l'/2'/3'(l-/l)(l-/2)(l-/3)] (3) 



P2P3 
P'l''2P'3 



where the scattering rate WI23 arises to fourth order 
in the bare two-particle interactions |31|. For low ener- 
gies and spinless fermions, Pauli principle ensures that it 
takes the universal form 



'P2){Pl ~P3.){P2~Pz) 



{p'i-pWi -pW2-P3))'Sp.,Pfm - Ef) (4) 

where Pi(/) = Pi + P2 + Pa and E^j^ = £1 + 62 + 63 
are the total momentum and energy of the three scatter- 
ing particles before (after) the collision, respectively. A 
simple dimensional analysis allows to identify a charac- 
teristic length scale of equilibration at the QCP (/i = 0) 
by setting typical momenta to \/2mT 



t-eq 



2Wm^L'^ 
(27r?i)4 



(2mT) 



13/2 



(5) 



Measuring all length scales in units of icq and all mo- 
menta in units of \/2mT allows to scale out the parame- 
ters W , m and T and the only remaining parameters are 
L/^eq, eV/T and p/T. We have checked both numerically 
and analytically that close to the QCP Hartree-Fock po- 
tentials (not included in Eq. ([!])) can be neglected. 

For a numerical solution of the Boltzmann equation ([T]) 
it is important to avoid discretization errors leading to a 



violation of conservation laws. We therefore use a con- 
servative splitting method following Ref. [1], see supple- 
ment [3^ , to solve the time-dependent Boltzmann equa- 
tion until a steady state has been reached. For the linear 
response calculation we use a linearized collision integral. 

Conservation laws. — Three conservation laws gov- 
ern transport in long quantum wires: charge, energy 
and momentum conservation. The corresponding cur- 
rents are the charge current, jc = ej^p'^pfpj the energy 
current, Je ~ J2p^p^pfp^ and the momentum current, 
jp w J2pP'^pfp- latter can be identified with pres- 

sure. For sufRciently long quantum wires and far away 
from the contacts the system will reach locally equilib- 
rium with the distribution function 



(6) 



parametrized by three space-dependent Lagrange param- 
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FIG. 1: (Color online) Electron density, n{x) = no{l + Sn{x)), 
at the QCP (n — 0) calculated from a solution of the Boltz- 
mann equation using (a) a linearized collision integral for 
L — lOO^cq and V ^ and (b) the non-linearized collision 
lO^oq and larger voltages, eV/T — 0.4 and 




integral for L 

eV/T = -0.6 (inset: eV/T = 0.4) 
T{x) - T and u{x) (in units of T and ^277 



The inset shows i-i{x) 
obtained from 



fitting the local charge, energy and momentum densities to 
Eq. ((6]). While in the linear response regime, L <g \iv\, there 
is a linear voltage drop across the wire, the voltage (and the 
density) drops predominantly close to one of the two contacts 
for L 3> \£v\- Note that due to a finite drift u the chemical 
potentials close to the contacts do not match the chemical 
potential /i ± eV/2 in the leads (shown as separate dots at 
X — +I//2 in the insets). 
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eters fJ,{x), T{x) and the velocity u{x) reflecting the three 
conservation laws. For the distribution function ^ one 
can calculate the corresponding equilibrium currents j^'i, 
j^^ and as function of /i, T and u. 

Voltage drop. — Fig.[I]shows the density profile and the 
local chemical potential (insets) of long quantum wires, 
L 3> ^cq, obtained from our Boltzmann simulations (to 
a good approximation, ^{x) can be measured by tunnel- 
ing contacts (l^l)- For small V (Fig. [1^) there is both a 
linear drop of the chemical potential along the wire and 
a finite jump directly at the two contacts (the separate 
points at a; = ±L/2 show ^ in the leads) on a length scale 
set by the geometry of the contact (set to zero within 
our approach). This jump is also present for larger V 
(Fig. [T]d), where, however, the linear voltage drop is ab- 
sent. Surprisingly, there is instead a large asymmetric 
voltage drop which occurs only close to one of the two 
contacts. This behavior occurs for sufficiently long wires 
not only directly at the QCP but also away from it. 

The qualitative difference between small and larger 
voltage can be understood from a simple argument based 
on matching currents. The steady state for F > is 
characterized by the three currents jc, js and jp. From 
the three equations ja = Ja^, a = c,E,p, one can, for 
sufficiently long wires, determine the three parameters 
/i, T and u which will be constant along the wire as 
ja — const.. Therefore, for a sufficiently long wire and 
finite V, a voltage drop can occur only close to the con- 
tacts. In the linear response regime, i.e. for small V , the 
situation is, however, different. By setting only one of 
the three parameters, u, to zero, two currents, j^'^ and 
j^^, vanish in equilibrium. As both and are linear 
in u, their ratio jc^/j'^ is - in the limit of small V - 
fixed by the average and T. This is used below when 
calculating the linear-response conductance analytically. 

To develop an approximate analytical theory valid in 
both regimes, we consider small, but finite voltages V 
and parametrize fx.p by 



fx.p — fx% ~^ Sfx,p, 



(7) 



where dfx,p accounts for deviations from local equilib- 
rium. Here is convenient to determine fJ.{x), T{x) and 
u{x) from the two equations jc — jc"^ and jp — j'^ while 
the third parameter is fixed by fitting the local density 
nix) = 71°'^ (a;). 

By linearizing the Boltzmann equation ([T]) in (5/, one 
obtains that 5f is proportional to dn/dx. For the energy 
current, one therefore obtains 



dx 



JE = j^E' + — const. 



(8) 



where _D is a certain thermoelectric diffusion constant 
which for small voltages is approximately constant across 
the wire. Using Galilei invariance which implies u = 
jl'^/en"'^ = jc/en, we obtain j"^ = 



V the last term can be neglected and one can linearize 
the density n = no + Sn to obtain 



•^jpjc 



Sn 



^d6n 
ox 



' JE 



'^jpjc 

2enQ 



~ const. 



(9) 



(10) 



This equation introduces a new length scale 

which diverges for V ^ as jc vanishes in this limit 
while D, no and jp remain finite. For \eV\ <C T and 
/i = 0, i.e., at the QCP, a simple dimensional analysis 
gives ^ ^cq-^- 

For L < \£v\, one obtains from Eq. ([9]) ^ = const. 
and therefore a linear drop in density and local chemical 
potential as in our numerical results, Fig.[T^. In the other 
limit, L 3> \lv\, Sn obtains an exponential x dependence 

-i/2" 



n{x) = nL + (nii - n^) exp 



(11) 



with nj^/ji « n(=pL/2). The direction of the current and 
the sign of D determines whether the drop of density 
and voltage occur at the right {ly > 0) or left (£y < 0) 
lead, see Fig. [T|d. This shows that the strongly asym- 
metric drop of voltage arises from a thermoelectric effect 
captured by the simple hydrodynamic equation Q. 

Linear response regime. — Interestingly, it is possible 
to calculate in the linear response regime the quantum 
critical conductance for long quantum wires {£oq ^ L ^ 
\£v\) analytically. We use the approach developed in 
Ref. [H, H (where only T < fi was considered) and keep 
track of the change of the charge and energy current 
carried by right-moving electrons with p > 0, jc {x) = 
eT,p>o^kfx.,p and j§{x) = Y.p>o^P^pfx,p, respectively. 
We use that far away from the contacts the distribution 
function obtains local equilibrium, f^^p ~ Z"'^, described 
by Eq. ^ with T{x) = T + ST{x), n{x) = fJ. + Sfi{x) 
and u{x). Current conservation implies u{x) — const. 
Furthermore, the ratio ri = js/jc can easily be calcu- 
lated as it is independent of u, i.e. just a function of 
^-independent parameters. The condition of constant 
momentum current fixes within linear response another 
ratio r2 = dxST / dxSfjL. This allows to calculate the ratio 

= dxj§/dxj^ as the remaining unknown dxSfj, drops 
from this expression. To leading order, ri,r2 and r^ are 
space independent functions of /i and T (calculated in 
the supplementary material ^3Qj ) . Finally, one identi- 
fies [ij the difference in the charge (energy) current of 
the interacting- and non-interaction system as the total 
change in the right-moving charge (energy) current along 
the wire, jc = j° + J dxj^dx {js = j% + f dxj§dx), 
respectively. If we now assume (as we checked numeri- 
cally), that these integrals are dominated by their bulk 
contribution, we obtain the equation 



"■'"^ . For small 



JE 



Jc - fc 



(12) 
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from which one can calculate directly jc- Combining all 
results |30| . we find for the linear-response conductance 
up to corrections of 0(4q/i), i-e. for £cq<S^ L <^ \£v\ 



G{z) 



ao{z)a2{z) - al{z) 



h a2{z) + aQ{z)K'^{z) — 2q!i(z)k(z) 



(13) 



where z = ^/T, (...) 

and ak = (e)., ^ = ^J^f- 
/i = 0, this gives 

Gqcp _ 



nde(...)^,with/o 



^ - 21n^2 

6 



e'^/h I 9 C"(3/2) I 6 C(3/2) 1 o 

' 3 + 8(Hl/2) + y2C(l/2) ^'^^ 



l+e« 

At the QCP, i.e., for 



0.420 (14) 



with C(a;) the Riemann zeta function. Gqcp is about 
16% below the non-interacting result e^/2/i. 

Fig. [2] displays the linear response conductance as 
functions of pi/T for non-interacting (see below) and 
fully equilibrated electrons which have a clearly differ- 
ent shape. Our analytical formula ([T^ fits very well the 
numerical result (symbols). 

Shorter wires. — Upon lowering T, i^q rapidly in- 
creases, see Eq. ([S])- For quantum wires, where L/i^q ^ 
1, one can neglect the effects of equilibrating interactions. 
Half of the voltage drops at the left and right contact, re- 
spectively, and there is no voltage drop inside the wire as 
fx,p = fp is independent of x. For jc one obtains the well- 
known non-interacting result j° — ^ In }i}]g(^-ev/2)/T 
The conductance plateau transition in linear response 
is therefore described by G{^l/T) = Go/(l + e"'^/'^) 
while at fi — the current is for arbitrary eV/T given 
by Jo — GoV/2. In the inset of Fig. [2] we have cal- 
culated numerically the crossover from the interacting 
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FIG. 2: (Color online) Conductance of fully equilibrated 
(solid line) and non-interacting electrons (dashed line) in the 
linear response regime (^eq <Si L <^ l^v|)- Numerical re- 
sults (symbols) agree with Eq. (|13p . Inset: Upon lowering 
T, i'cq grows rapidly and a crossover from the equilibrated 
to the non-interacting conductance is observed for ^oq ~ L 
{L= 104q,At = 0). 
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FIG. 3: (Color online) Non-linear conductance, jc/V obtained 
for L — 7.5£cq from a numerical solution of the Boltzmann 
equation. For eV/T ^ a linearized collision integral was 
used. Inset: Within our numerical precision, there is no finite 
size dependence of the non-linear conductance for L ^ ^eq {L 
measured in units Icn )■ 



quantum critical conductance (|14|) to the non-interacting 
one, which occurs upon lowering T when ^cq ^ L. 

Nonlinear response. — Fig. [3] show the nonlinear con- 
ductance jc/V at the QCP, i.e. for /i = 0. It inter- 
polates between the linear-response value (fH)) and the 
non- interacting result obtained for \eV\/T — c». For 
\eV\ > T and ^ = all states withp > and ep < \eV\/2 
are occupied. As this is also an equilibrium distribution 
function with u = \J\eV\l Am and [i = mv? /2, collisions 
have no effects in this limit. 

For small V the nonlinear conductance appears to be 
non-analytic. 



Jc{V) = GqcfV + -f\V\V + . 



(15) 



which can be traced back to the asymmetric voltage drop 
for L ^ As riR — til in Eq. (jlip varies linear in 
the density in the center, n(0) ~ max(?ifl,ni) according 
to Eq. p4|) . obtains for L ^ \lv\ a correction propor- 
tional to As jc ~ eun, this implies a correction 
proportional to |I^|I^ to the current as soon as L » 
Due to numerical problems, we were not able to obtain 
reliable numerical results in the small- 1^ regime L \iy\ 
where we expect a rounding of the non-analytic correc- 
tion. Overall, the finite size corrections to the non-linear 
conductance are smaller than our numerical resolution 
for L:^ icq, see inset of Fig. [3l 

Outlook. — While our results have been derived only 
for spinless fermions with short ranged interactions, we 
expect that our main qualitative results are also of di- 
rect relevance for quantum wires made of electrons with 
spin and long-ranged Coulomb interactions. For these 
systems, the hydrodynamic equation ^ should also be 
valid (with strongly modified parameters). The highly 
asymmetric voltage drop predicted by us will probably 



5 



be even much easier to observe, as the stronger interac- 
tions imply that the regimes L » ^eq and L ^ \£v\ are 
much easier to reach. 
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Irrelevance of short ranged interactions at QCP 

The imaginary-time action describing spinless fermion 
with short ranged interaction at the QCP (/i = 0) takes 
the form 

(S16) 

The derivatives in the interaction term reflect Pauli's 
principle, ip'^ — 0. From this action one can directly read 
off the scaling dimensions of the field, [-0] ~ L~^/'^, and 
the (imaginary) time, [r] ~ L^, for [x] ^ L. This sim- 
ple power-counting shows that the interaction strength 
scales to smaller values, V — > V/X, as the distance of 
Fermions is increased, x — > Xx, or as temperature is 
lowered, T T/X^. Therefore the interactions are ir- 
relevant at the QCP. This result also implies that the 
quasiparticle picture and therefore the Boltzmann equa- 
tion becomes asymptotically exact upon approaching the 
QCP. 



and 6fx,p is linear in V , and linearize the collision inte- 
gral in 6f. For the numerically more demanding calcula- 
tions employing the full collision integral we used meshes 
with 22 and 42 points in momentum space and various 
different homogeneous and inhomogeneous discretization 
of space with ^100 grid-points. We checked that our 
results are independent of discretization. 



Linear conductance in long wires 

A method to calculate the linear conductance in long 
wires icq <C L <C |^y|, based on conservation laws and the 
distribution of fully equilibrated electrons Eq. (6) in the 



(6) ii 

main text, was originally introduced in Refs. [l|, y. As 
described in the main text, we need to calculate the ratios 
ri = jE/jc, r2 = dxST/dxS^j, and = dxj§/d^j^. The 
current jc in response to the applied voltage (to linear 
order V) is then found from 



Jc-fc 



(S17) 



Conservative splitting method 

For the numerical implementation of the Boltzmann 
equation we follow Rcf. 3] and solve the time-dependent 
Boltzmann equation, by splitting time evolution into a 
free flow and a relaxation stage. The advantage of the 
splitting procedure is that the distribution obtained after 
a relaxation step can be corrected such that conservation 
laws are fulfilled in each collision, see Ref. Q for details. 

Free fiow and relaxation steps were implemented by 
a finite difference scheme, i.e. by discretizing two- 
dimensional phase-space. We used a first order implicit- 
explicit upwind scheme to model the free propagation 
step, and an implicit scheme for the relaxation step Q. 
The steady state is reached when currents of the con- 
served quantities, jc, jp and Je are constant along the 
wire. For the linear response calculation we parametrize 
.fx,p — fp+Sfx.p, where is the distribution at zero bias 



where jj? and are the (linear response) charge and 
energy currents of non-interacting electrons which is di- 
rectly described by Eq. (2) of the main text, 



•0 e'V 



Je 



^ (/iao + Tai) , 



(S18) 



Here and in the following 
and z = fi/T. 



With these definitions we first calculate ri. Using the 
equilibrium distribution Eq. (6) of the main text we find 



n — ^ — - (fi + Tk) , K 

Jc e 



(S19) 



To calculate r2 we use that momentum conservation 
implies homogeneity of momentum-current in the steady 
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state. The latter can again be calculated with help of 
Eq. (6) given in the main text by expanding T(x) = 
T + ST{x), ^{x) = ^ + Sfi{x) form small V. To linear 
order in V, i.e., for small ST and we find 

const. — jp{x) — nSfi{x) + nKST{x) + const. ^ (S20) 

resulting in 

dJT 



r2 



dxS^ 



(S21) 



To obtain we can directly use the definition of 
and j§ given in the main text combined with the equi- 
librium distribution function (6) of the main text 



OxJc = -[-OxOfJ- + ——dxOT 
n h 



(S22) 



and 



= + ^ {al^J^l + a2dJT) . (S23) 



For the ratio r2 we therefore obtain 

^ _ Q^Je _ f-^ ^ T aiK - Q!2 
dx.i^ e e aoK- ai 



(S24) 



Inserting above expressions into (|S17p and solving for jc 
gives Eq. (13) of the main text. 



Measuring voltage profiles 

To measure the voltage drop across a quantum wire 
one can, for example, use a weakly coupled tunneling 
contact realized by the tip of a scanning tunneling micro- 
scope. Assuming a constant tunneling matrix element M 
and a constant density of states uq of the tunneling tip, 
the charge current Ic{x) and the energy current Ie{x) 
through the tip located at position x are given by 

Ic{x) = ^\M\' [dk{fx^,~ff%{e,)) (S25) 



Ie{x) 



Here fx,p is the distribution function of the wire and 
/*'P(e) is the Fermi distribution describing the occupa- 
tion of the states in the tip. The latter is parametrized 
by the chemical potential fl and the temperature T. 

The local chemical potential fJ-{x) of the quantum wire 
and the local temperature T(x) of the wire are now ob- 
tained from the condition that particle and energy cur- 
rents flow only if there is a difference in the chemical 
potential and temperature. T{x) and fJ.{x) are therefore 
obtained from the condition 

T{x)=f, tJ.ix)=)l for = /^(x) = (S26) 



^ 0-2 
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FIG. S4: Local chemical potential fJ^{x) (blue) and temper- 
ature T{x) = T + 5T(x) (red) (in units \eV\) for a wire of 
length L = lO^cq and voltage eV /T = 0.4. Symbols: Profiles 
obtained using Eq. (|S26|) . i.e. for voltage contacts in the tun- 
neling regime. Solid lines: Corresponding profiles obtained 
from fitting local charge, momentum and energy densities in- 
side the wire to an equilibrium Fermi distribution with de- 
pends not only on fj, and T but also on a local velocity u. 

Note that this definition can be used for distribution 
functions far from equilibrium (and the usual result is 
obtained in equilibrium). Eq. (jS25P implies that the 
local chemical potential and temperature can directly be 
obtained from the local charge- and energy density of the 
system. 

In Fig. lS4l we show chemical potential and temperature 
profiles obtained from the definition (jS26|) for a wire of 
length L = 104q with eV/T = 0.4 at the QCP {n = 
0). The separate points at ±i/2 denote the chemical 
potential in the two leads. 

In the insets of Fig. 1 of the main text we also show 
/i(x) and T{x) but in this case we use a different defi- 
nition of these quantities, which is more appropriate to 
illustrate our analytical arguments. In the main text, we 
fit the local charge, momentum and energy densities to 
the equilibrium distribution given by Eq. (6) of the main 
text which depends not only on /i and T but also on the 
average velocity u. In contrast, according to Eqs. (|S25p 
and (|S26|) . the local charge and energy density are fitted 
to a distribution function with w = as momentum is 
not conserved in the tunneling tip. 

For the range of applied voltages discussed in this pa- 
per, both methods lead to nearly identical profiles. This 
is shown in Fig. IS41 where temperature and chemical- 
potential profiles obtained from Eq. (|S26I) (symbols) are 
compared to the corresponding curves from the main text 
(lines). 
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